Quasi-stationary trajectories of the HMF model: a topological perspective 
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We employ a topological approach to investigate the nature of quasi-stationary states of the Mean 
Field XY Hamiltonian model that arise when the system is initially prepared in a fully magnetized 
configuration. By means of numerical simulations and analytical considerations, we show that, along 
the quasi-stationary trajectories, the system evolves in a manifold of critical points of the potential 
energy function. Although these critical points are maxima, the large number of directions with 
marginal stability may be responsible for the slow relaxation dynamics and the trapping of the 
system in such trajectories. 
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The so called Mean Field XY Hamiltonian Model 
(HMF model) Q has received great attention during 
the last years in the statistical mechanics community, 
mainly because of the richness of its dynamical behavior 
@ 0, II 0, H 0, H- The model is defined by a set of A^ 
particles, or rotators, moving on a unitary circle. The 
dynamics of the system is ruled by the following Hamil- 
tonian: 
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Here 6i represents the rotation angle of the z-th parti- 
cle (with 9i 6 (0, 2n]) and Pi its conjugate momentum. 
This model can be considered as a kinetic version of the 
mean field XY magnetic model with ferromagnetic in- 
teractions. From a thermodynamical point of view the 
model is extremely simple, in contrast to its rich and 
still not well understood dynamical behavior. On one 
side, being a mean field model, its equilibrium thermo- 
dynamics can be exactly solved in the canonical ensem- 
ble, yielding a second order ferromagnetic phase transi- 
tion. On the other side its relaxational dynamics is very 
complex and equilibrium is not easily attained from an 
important set of initial conditions 0,0,111]. The origin 
of this kinetic complexity is elusive and some similarity 
with the phenomenology of disordered systems has been 
advocated Nevertheless, the existence of this re- 

lation is not at all obvious. In first place, there is no 
imposed disorder on the Hamiltonian. Second, the cou- 
plings are all ferromagnetic, avoiding in this case any kind 
of structural frustration. Finally, the infinite range of the 



interactions further simplifies both the dynamics and the 
thermodynamics of the model, avoiding any topological 
consideration on the structure of the lattice where par- 
ticles are located. On the other hand, and also due to 
the infinite range of the interactions, its dynamics can 
be efficiently integrated in computer simulations. There- 
fore, this is an excellent prototype for analyzing the mi- 
croscopic dynamics of a finite system close to a critical 
point. 

The microcanonical simulations of the HMF reproduce 
many of the anomalous critical behaviors observed in nu- 
clear and cluster fragmentation processes. In particular, 
depending on the initial preparation of the system, it is 
possible to verify the existence of negative specific heat 
curves, qualitatively similar to those observed in recent 
fragmentation experiments in small clusters (see for in- 
stance 0| and references therein). But the interest on 
this model largely exceeds this motivation. The exis- 
tence of Quasi-stationary solutions whose life-times di- 
verge |2,l3|, in the thermodynamical limit N — > oo, has 
raised the question on whether it is possible or not to 
construct a measure theory able to predict the station- 
ary values of physical observables in long standing out 
of equilibrium regimes [ijj- Furthermore, the existence 
of a glassy-like relaxation dynamics 0, 0] along quasi- 
stationary trajectories, has opened new challenging ques- 
tions on the origin of such unexpected behavior for an un- 
frustrasted non disordered model. Finally, in the last few 
years, the HMF model has been used, as a paradigmatic 
example , fo r the study of the so called Topological Hy- 
pothesis |l2 ; LUiLLJ] which asserts that phase transitions, 
even in a finite system, can be identified by searching for 
drastic topological changes in the submanifolds of the in- 
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teraction potential. 

Concerning the HMF thermodynamics pj, the usual 
microcanonical and canonical calculations predict that 
the system suffers a second order phase transition at T c — 
1/2 (which corresponds to an internal energy per particle 
U c /N — 3/4), from a low temperature ordered phase to 
a high temperature disordered one. Actually, one can 
associate to each rotator a local magnetization: m, = 
(cos0j,sin0j), and then define the order parameter of the 
transition as the global magnetization: 



M = (M x ,M y ) = ±J2'' 
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At the critical energy, M = \M\ vanishes continuously as 
the system is heated from the ordered phase. 

Concerning the dynamics, when the system is prepared 
very far from equilibrium, for energies just below the crit- 
ical one, the system gets trapped into quasi-equilibrium 
trajectories. These trajectories are characterized by time 
averages of one-time observables which reach, after a 
rapid initial transient, almost constant values which do 
not coincide with those predicted by microcanonical or 
canonical ensemble calculations Q. The average time 
that a system of size N remains in a quasi-stationary 
trajectory grows with N Therefore, if the system 

were infinite, it would remain there forever, without ever 
reaching true equilibrium. An even more surprising sce- 
nario appears when one considers the relaxation of the 
two-time correlation function C(t, t') (either in the whole 
hase space j3| or considering only the momenta space 
). The explicit dependence of C on both times t and t' 
indicates the loss of time-translational invariance, proper 
of equilibrium states, and the appearance of memory ef- 
fects, a phenomenon usually called aging. The scaling law 
of the two-time autocorrelation functions 0, 0] is quali- 
tatively similar to that observed in some real spin glasses 
[lfij . Nevertheless we will show that the physical mech- 
anisms behind the quasi-stationary states of the HMF 
are completely different from those present in disordered 
systems. 

In this work we will show, through numerical simula- 
tions, that the complex nonequilibrium quasi-stationary 
regime observed just below the critical energy can be in- 
terpreted from a topological point of view. Our analysis 
will focus on the topology of the surface defined by the 
potential energy in the configuration space. The HMF 
potential energy, as well as any Curie- Weiss like poten- 
tial, can be written in terms of the order parameter of 
the system: 
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Note that the potential energy per particle V/N takes 
values in the interval < V/N < v c = 1/2. The lower 



limit corresponds to the case of the fully ordered config- 
urations (hence M = 1) and the upper bound to a com- 
pletely disordered configuration. The configuration space 
manifold M is an A^-dimensional torus parametrized by 
the N angles 6i, The critical points (CPs) of Ad are 
those points for which all the N derivatives of V/N van- 
ish, i.e., d(V/N)/d9i = 0, for i = 1, . . . ,N, Making use 
of the infinite range of the interactions, one can write the 
derivatives of the potential in terms of the two compo- 
nents of the order parameter, namely, 
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Furthermore, CPs can be classified according to the 
eigenvalues of the Hessian of V/N, that for the HMF can 
be written as H = D + B where 
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In our work we use the following protocol: starting 
from a far-from-equilibrium configuration, we integrate 
numerically the set of Hamilton equations of the system 
using a fourth order symplectic method with a very small 
time step (typically dt = 0.01). Along the trajectories we 
evaluate, at each time step, the modulus of the N deriva- 
tives of V/N given by <@J and identify the maximum over 
i = 1, . . . , N, through 

d{V/N) 
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Then, each time that A = 0, it means that the system 
reaches a CP (actually, due to the finiteness of time step, 
A is never exactly zero, but it gets closer as dt decreases). 

Let us first analyze the behavior of the system in the 
disordered phase. It is important to stress that, although 
the system ultimately relaxes to equilibrium, a slow re- 
laxation has been observed when starting very far from 
equilibrium This is probably due to the fact that 
rotators move almost freely and trajectories are weakly 
chaotic In Fig. we plot V/N and A as a function 

of t, for U/N = 10, in the disordered phase well above 
the critical energy U c /N — 3/4. The system has been 
initially prepared in a "water-bag" configuration, with all 
the rotators aligned along the x axis (9, = 0, for all i) 
and the momenta drawn from a uniform distribution (ac- 
tually we used regularly spaced momenta 0])- Fig. 
indicates that the system periodically visits CPs of the 
potential, corresponding to V/N =1/2 (hence M = 0). 
In fact the observ ed p eriod is of the order of the mean 
period of rotation [16] . 

Let us now compare these results with those obtained 
in the low energy phase, just below the phase transi- 
tion. In this case, quasi-stationary solutions emerge, dis- 
playing a kind of glassy-like dynamics characterized by 
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FIG. 1: Time evolution of (a) A, the largest modulus of 
the derivatives of V, and (b) V/N, the potential energy per 
particle, for a system of 500 particles and U/N — 10, well 
above the phase transition. The system was initially prepared 
in a regular water-bag configuration. The dotted line in (b) 
corresponds to the equilibrium value in the thermodynamic 
limit. 



energy phase, the time intervals At elapsed between two 
successive CPs do not present any pattern of periodic- 
ity. On the contrary, the system visits the CPs in an 
apparently disordered way. The probability distribution 
function (PDF) -P(At) of time intervals between CPs is 
shown in Fig. 03 P(At) can be reasonably fitted by a 
power law decay. Other striking features of the dynam- 
ical behavior can be noted in Fig. |2 First, the system 
is initialized in a configuration at the bottom of the po- 
tential energy V/N = and it goes almost abruptly to a 
region with a mean potential energy per particle larger 
than the equilibrium mean potential energy, and stays 
around this level during the whole time span of the sim- 
ulation. In fact it goes close to the top of the potential 
energy landscape, V/N » 1/2, tries to escape downhill 
but uses the kinetic energy gain to attain again CPs at 
the top. Somehow the system is not able to relax from 
this level to the equilibrium level during very large time 
scales. The configurations sampled correspond to the 
quasi-stationary states and the system stays there dur- 
ing time spans which scale with the size N, being trapped 
forever in the thermodynamic limit. One can say that, 
in configuration space, quasi-stationary states are always 
near CPs of the landscape. 
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FIG. 2: Time evolution of A (a) and V/N (b) for a system of 
500 particles and U/N = 0.69, just below the phase transi- 
tion. The system was initially prepared in a regular water-bag 
configuration. The dotted and dashed lines in (b) correspond 
to equilibrium and metaequilibrium values in the thermody- 
namic limit, respectively. 



weak chaos, non-Gaussian velocity distributions and sub- 
aging, as mentioned above. In Fig. |2j we plot V/N and 
A vs. t for U/N — 0.69 and water-bag initial conditions. 
Here we verify that again the system sequentially visits 
one CP of the potential energy after another, also corre- 
sponding to M = 0. However, at variance with the high 
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FIG. 3: Probability distribution function P(Ar) of the time 
intervals between two consecutive critical points, for the same 
system of Fig. [2] 



We next address the effect of system size. The same 
qualitative behavior can be observed, with a bigger nu- 
merical effort, for larger systems. The departure from 
the CPs (measured by the height of A between consecu- 
tive CPs) decreases as N increases, but the distribution 
of Ar remains unaltered. From Eqs. (0) and Q, it is 
straightforward to see that at the critical level v c = 1/2 
the CPs are continuously degenerate. Consequently as N 
grows and fluctuations outside the critical level diminish 
the system wanders more and more inside this manifold 
of CPs. As the energy decreases, the CPs are visited more 
sparsely, down to U/N ~ 0.67, where the system stops 
wandering among CPs. It is noteworthy that spatially 
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homogeneous states lose Vlasov stability approximately 
at that energy (l7j . 
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FIG. 4: Time evolution of A and V/N for a system of 500 
particles and U/N — 0.69, just below the phase transition. 
The system was initially prepared in a close to equilibrium 
configuration. The dotted line corresponds to the equilibrium 
value of V/N. 



A completely different scenario emerges when the sys- 
tem is prepared, at a given energy, in an almost equi- 
librated initial configuration. The plot of A vs. t pre- 
sented in Fig. 31 clearly indicates that, in this case, the 
system does not wander among CPs. Its potential en- 
ergy per particle (and then also its temperature) rapidly 
starts fluctuating around their canonical values, showing 
strong finite size effects. 

A crucial piece of information comes from the stabil- 
ity of the visited CPs. Since capturing the exact time 
at which the dynamics passes through a CP is a difficult 
task, numerical evaluation of the Hessian at CPs arising 
from the dynamics may lead to wrong estimates. There- 
fore it is important to perform some analytical calcula- 
tions. In order to do so, let us recall that, both in the high 
energy phase and in the quasi-stationary states, the CPs 
correspond to V/N =1/2 (see Figs. 1 and 2), hence they 
are points of M. with zero magnetization. Moreover, the 
distribution of angles at those points is approximately 
uniform 0|. A configuration with these characteristics 
for which analytical calculations are possible, consists in 
regularly distributed angles in the interval (Q,2tt], i.e., 
6 k = 2nk/N, k = 1, N for even N. If M = 0, from 
Eq. (|3|l. the Hessian of the potential energy is H = B. 
Then, for the regular configuration we have 



H = B= — (1 + A), (7) 

where Am — cos(27r(fc — l)/N). This circulant matrix can 
be diagonalized in Fourier space, yielding the following 



eigenvalues of the Hessian matrix 



Hi = -jpl cos(7i7) + 2 ^2 cos(2irrl/N) cos(2irr/N) j , 
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for 1 < I < N. Thus, we obtain B x = H N -i = -l/(2N) 
while the remaining eigenvalues are null. This means 
that at these CPs there are two unstable directions and 
N — 2 marginal ones. This picture remains valid for more 
general situations than those restricted to the particular 
regular cases. In fact, we have verified that when angles 
are randomly chosen in the interval (0, 2ir] there are two 
eigenvalues with values close to —1/(2N) and the remain- 
ing A^ — 2 vanish. Finally, the eigenvalues calculated from 
the configurations of the dynamics very close to CPs are 
also consistent with this picture. >From this analysis 
it is clear that in the high energy phase and also in the 
quasi-stationary low energy states, the systems wanders 
in an almost flat landscape. 

In the high energy phase, energy is mainly kinetic, 
leading to ballistic behavior in the flat landscape of 
CPs during long time scales. In the low energy, quasi- 
equilibrium regime, kinetic energy is comparable to the 
potential one. At this relatively low energies diffusion 
is slower than ballistic. Nevertheless, due to the flatness 
of configuration space superdiffusive behavior is observed 
1- 

Summarizing, we have seen, by means of numerical 
simulations, that the quasi-stationary trajectories ob- 
served in the HMF model can be interpreted in terms 
of the topological properties of the potential energy per 
particle V/N of the model. Starting from a far from 
equilibrium configuration with V/N = 0, the system ini- 
tially decreases as much as possible its kinetic energy and 
settles at the flat top of its potential energy. The sim- 
ulations confirm that along the quasi stationary states 
the system wanders among different CPs with only two 
negative directions and N — 2 marginal ones. Moreover, 
CPs at the upper critical level are continuously degen- 
erate. This suggests that once inside this critical sub- 
manifold the system cannot go out easily and relax to 
equilibrium. A probable scenario in the thermodynamic 
limit is that, as the fluctuations of the potential energy 
close to the critical level v c — 1/2 go to zero, the system 
keeps wandering continuously inside the manifold of CPs 
and consequently remains forever out of thermodynamic 
equilibrium. 
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